Cosmic PeV Neutrinos and the Sources of Ultrahigh Energy Protons 
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The IceCube experiment recently announced the first detection of multi-PeV neutrinos, which cannot easily 
be explained by atmospheric or cosmogenic fluxes. We examine whether these neutrinos originate from the same 
sources as ultrahigh-energy cosmic rays. We find that producing the requisite neutrino flux through photopion 
production in the source leads to a proton flux at the level of cosmic-ray data at ~ 10 18 eV in a constrained 
scenario where neutrinos only arise from ir + decays. In more general cases the proton yield is much lower, 
requiring a dominant class of accelerator that allows cosmic rays to escape without significant losses. 
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Introduction. — High-energy astrophysical neutrinos have 
much to tell us about the most extreme environments in the 
Universe; however, finding them is a difficult endeavor 03- 
3|. Colossal detectors are required, such as IceCube [4|, that 
can observe the tracks of muons produced in charged-current 
Vft/Vfj, scattering or showers (cascades) induced by various 
channels (as we discuss later). The recent observation of two 
PeV-energy shower events by IceCube may represent the dis- 
covery of such neutrinos [5|, as atmospheric PeV neutrino 
fluxes are low ||6l|7). The most likely astrophysical mecha- 
nism is photopion production by protons on a photon back- 
ground, p 7 — > N it, leading to neutrinos via pion decays. 

An example of this process is the suppression of ultrahigh- 
energy cosmic-ray (UHECR) proton fluxes at > 10 19 5 eV due 
to the cosmic microwave background (CMB), the GZK ef- 
fect |8] [9j. The measured UHECR spectrum displays a down- 
turn near this energy lfl0l - [T3l . However, the > 10 18 eV neutri- 
nos resulting from the decays of pions produced lfl4l4T7ll are 
too energetic to explain the IceCube events, while the flux of 
lower-energy neutrinos from decays of neutrons is too small. 

The cosmic infrared/optical background allows for photo- 
production by lower-energy protons in propagation, leading 
to lower-energy neutrinos [ 18-20 1. However, Fermi measure- 
ments of gamma-ray absorption in blazar spectra now indi- 
cate a low level of the < 10 eV diffuse photons ETTl needed to 
photoproduce with < 10 17 eV protons to yield ~ 10 15 eV pi- 
onic neutrinos. Moreover, unless magnetic fields lead to suf- 
ficient synchrotron losses in the subsequent cascades to sup- 
press gamma rays, a PeV neutrino flux cannot be produced at 
the required level that respects the isotropic gamma-ray back- 
ground [ 22 1 . This suggests that, if the neutrino flux is cosmic, 
it likely arises from within some class of sources. 

Our goal is to determine whether the IceCube PeV neutri- 
nos share a common origin with UHECRs in the ~ 10 18 eV 
range where the composition is inferred to be light Il23lj25l . 
The UHECR protons observed at Earth must have been able 
to retain sufficient energy upon leaving their acceleration sites 
11261 . However, the magnetic fields required to contain the par- 
ticles during the acceleration process can lead to severe adia- 



batic losses if they cannot be promptly escaped. 

If this is to be achieved through photohadronic interactions 
producing neutrons that later decay to protons outside the 
source, then there must be an accompanying flux of pionic 
neutrinos (for this reason, IceCube limits disfavor gamma-ray 
bursts [27 -30 1). We use the neutrino flux implied by the Ice- 
Cube events (for an assumed source spectrum, evolution with 
redshift, flavor ratios, and neutrino oscillations) to normalize 
the outgoing proton spectrum to determine whether this is the 
mechanism operating in ultrahigh-energy sources. 

Neutrinos and IceCube. — We begin by constructing a 
source neutrino spectrum that results in fluxes at Earth that 
yield shower rates that peak at the IceCube energy range us- 
ing a smoothly-broken power law (with E in GeV) as 
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FIG. 1: Total neutrino fluxes (normalized to yield two > 1 PeV show- 
ers in IceCube) in our models based on 7r + -only (which has no v\ top 
panel) and ir [i (bottom panel) decays, separated after oscillations 
into: v e , i> e (solid), v^, (dotted), and v T , u T (dashed). The Ice- 
Cube 40-string E~ 2 flux limit 1311 is also indicated. 
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where the slopes are a and /?, with r] = — 1 to give a smooth 
break at Ej,. We obtain the neutrino fluxes at Earth, cp ll (E u ), 
by integrating these spectra up to z max = 8 as 

c /•-»• dA^dEi W(z) ... 

V>v{E v ) = —\ ^73^37T^« Z : ( 2 ) 
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where dz/di = H (1 + z)[fi m (l + z) 3 + ^a] 1 ^ (with fi„ 
0.3, 51 a = 0.7, and H = 70 km/s/Mpc) and dE' v jdE v = (1 + 
z) accounts for redshift. We set the evolution in the rate of 
neutrino production, W(z) = 1, which conservatively bounds 
the required neutrino emissivity, and discuss alternatives later. 

Motivated by the model in ll32l for an E~ 2 accelerated pro- 
ton spectrum, we first use a neutrino spectrum with a = — 1, 
f3 = —2, and photopion opacity ~ 1 for Ei = 10 7 . We assign 
the initial flavor ratios resulting from equal numbers of 7r + 
and 7T~, as in sources with hard photon backgrounds. Neu- 
trinos result from 7r//i decays, ir + — > p^v^, p + — > e + v^v e 



and it" —> n~Pn, p~ — > e~v^v e , giving v e :Vn:v T = 1:2:0 and 
v e :v^:v T — 1:2:0. We account for neutrino oscillations using 
mixing parameters from 11331 . The lower panel of Fig.[T]shows 
the shapes of the arriving fluxes, broken down by flavor. 

We also consider the extreme in which only ir + are pro- 
duced, as with near-threshold photoproduction due to a soft 
background. Muons experience higher synchrotron losses 
than pions in the strong magnetic fields within the accelera- 
tion region of sources Il34ll35ll . We use this to obtain a model 
with no v component by assuming that muons are produced, 
but cool and decay only to low-energy neutrinos. The initial 
flavor ratios are approximated as those from the n + decay, 
Ve.Vp.VT =0:1:0 and v e \v^\v T =0:0:0. Here, we take pion 
cooling to result in a break at E\, — 10 7 , so that f3 = —3, as in 
the upper panel of Fig. [T] In both cases, we neglect the v e 
flux from neutron decay, which carries much less energy and 
peaks at energies lower by about two orders of magnitude. 

Rates of shower-like events in IceCube are principally de- 
termined by the normalized neutrino fluxes, effective vol- 
ume, and cross sections for neutrino-nucleon, <j v n and a^N, 
and antineutrino-electron, cr Pe , scattering. We use the total 
deep-inelastic scattering cross sections from 1 36 ] for charged- 
current (CC) and neutral-current (NC) scattering, and the av- 
erage inelasticity, (y(E v )), from 11371 . approximating (y) — 
0.25 at > PeV energies. 

The visible energy of the shower depends upon the inter- 
action channel. For NC events {v N —>v X), onc is identi- 
cal for all flavors, and (y) determines the fraction of energy 
imparted to a quark in the nucleon. The resulting hadronic 
shower results in less light than an equivalent-energy electro- 
magnetic shower by a factor that is a function of energy l38l 
that we assume to be /had — 0.9. Thus, working in terms 
of the electromagnetic-equivalent shower energy E em (as in 
Fig. [2j, with <7nc ~ 0.23, we have 

-^em,NC = fhadiu) E u — t/NC E v . (3) 

For CC v e events {v e N — > e X), we assume that the elec- 
tron deposits its entire energy, E e = (1 — y)E v , into an elec- 
tromagnetic shower. There is also an accompanying hadronic 
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FIG. 2: Neutrino event spectra in IceCube using the shower energy, 
E eln , from the fluxes in the top and bottom panels of Fig. [TJ Shown 
are: v e , z? e , v T , and v r charged current; all flavor (yx, vx) neutral 
current; and £/ e e channels yielding e and hadrons (included in v R 
line) and r. Normalizations are based on two total > PeV events. 



shower of energy (y)E v . We add these together to get a total 
effective visible energy per interaction as E cm . e = ( 1 — y ) E v + 
fhad{y)E u = q e E v , with q e w 0.95 (and similarly for u e ). 

The Glashow resonance, v e e — > W~— > X, is important 
near i? Pe w 6.3 PeV. W channels yielding quarks are purely 
hadronic (qo,q ~ 0.9). The e/r channels result in a neutrino 
carrying away most of the energy ((y) ss 0.25 [37|). 

Tau neutrino CC events ( v T N — > r X) typically have prop- 
erties intermediate between CC v e and NC events. One dif- 
ference is that a tau at PeV energies will travel ~ 50 m before 
decaying ll39l . with ~ 80% of decays 1331 involving channels 
that result in a shower, although with an escaping v T carrying 
away energy. Decays to electrons result in an electromagnetic 
cascade with an outgoing v e . Other decays involve multiple 
mesons, which result in a hadronic shower. Assuming all such 
decays to give hadronic -like cascades, including the initial 
cascade and assuming the two bangs to be indistinguishable 
for now gives E emt T = qr E v , with q T « 0.9. 

The spectrum of events for each channel can be given in 
terms of the electromagnetic-equivalent energy as fiOl |4D 
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where Na p is the molar density of ice. Using 27r sr, due to the 
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attenuation of upgoing PeV neutrinos, and V e g T ss 1 km 3 yr 
roughly matches the PeV IceCube exposure [5 |. 

We integrate the total vN, vN, and v e e (with the elec- 
tron density lower by 10/18) shower yields above 1 PeV, and 
equate to two events (assuming that all showers above this 
limit would be counted). For the n^fi^ flavor ratio model, 
we thus obtain the normalization for dN v jdE v using Eq. (JT|i 
as /o,7rp~4.9x 1CT 49 GeV _1 cm _3 s _1 , giving the total flux 
curve in Fig. [T] Our estimates using the da / dy distributions 
in 11371 agree at the ~ 10% level. The lower panel of Fig. [2] 
shows the resulting shower spectrum in the detector. 

Examining the importance of the Glashow resonance (see 
also Il42ll43l ). we note that D e e scattering results in the peaks 
at ~ 1.5 PeV and ^6 PeV, although the exact da/dy signifi- 
cantly broadens the former. Since the W width to hadrons is a 
factor of ~ 3 larger than to the sum of e/r, 6 PeV showers are 
the most likely and account for < 1 of the predicted events. 

In our £/-less scenario where only tt + are produced and 
muons cool to low energy prior to decaying, resonant D e e 
events do not occur. Using the flavor ratios from the ir + de- 
cay alone and starting again from Eq. ([T]i results in the shower 
signals in the upper panel of Fig. [2] To obtain two > 1 PeV 
events requires /o.tt = 1.4x 10 -48 GeV _1 cm" 3 s _1 . 

Muon neutrino CC scatterings (v^N — > /il) also produce 
a hadronic shower; however, a defining characteristic of such 
events is the resulting muon. Using the normalized neutrino 
spectra, and accounting for the ~ 20% of tau decays that re- 
sult in a muon, we find the expected number of > PeV muons 
with contained vertices to be ~ 1 in both cases, in accord with 
a lack of IceCube muon events in this range thus far. 

The cosmic-ray spectrum. — A long-standing hope is to 
determine the UHECR sources and ascertain the acceleration 
mechanism. By making a few simplifying assumptions, we 
relate our neutrino spectra to a proton flux with a normaliza- 
tion that is fixed by the IceCube data. 

tron with E n ~ 20 E v is produced corresponding to each tt^ 1 



We assume for the 7r ± u ± flavor ratio model that one neu- 



pair. Since six total neutrinos result from the n fi decays, 
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with neutron decay giving a source proton spectrum in the 
form of Eq. ^ with a = — 1, = -3, and E b w 10 83 due to 
neutron interactions in the source IT321 as assumed above. In- 
tegrating over 10 15 < E p < 10 21 

-3 



< E p < 10 21 eV, we find an emissivity 



3 - 8 -l:i x 10 44 erg Mpc" 



VI' 



(6) 



with uncertainties corresponding to the 68% Poisson confi- 
dence ranges [44 1 for two detected IceCube events. Introduc- 
ing a cutoff below 10 21 eV scarcely affects these energetics. 

The proton flux associated with the ir + scenario is taken to 
have the same spectral shape. However, (i) we must account 
for /o,7r ~ 3 /o,ir^; (m) fo.ir normalizes a flux from a single 
neutrino species, rather than the sum of six, so that the proton 
flux must be larger by an additional factor of six, thus giving 

£ P ,n ~ 6.7±|;1 x 10 45 erg Mpc" 3 yr" 1 . (7) 
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FIG. 3: The ultrahigh-energy cosmic-ray spectrum. Shown are the 
proton fluxes associated with our tt fx (dashed) and -n + (solid) 
neutrino flavor models, with normalizations obtained from the Ice- 
Cube neutrino observations, with 68% confidence bands based on 
two events. These are compared to HiRes-II 1 10], Auger 1121 . and 
Telescope Array 1 13 1 data. 



To find the expected proton fluxes at Earth from these spec- 
tra, three types of energy loss must be accounted for in propa- 
gation. At energies where photopion production on the CMB 
occurs (> 10 19 5 eV), 737 — > is the dominant loss chan- 
nel E). For E p > 10 18 eV, and below the CMB photo- 
pion threshold, resonant pair production on background pho- 
tons, p7 — > pe + e~, dominates |46|. This process has a large 
cross section, although each interaction removes only a small 
amount of energy. Finally, there is an energy-independent adi- 
abatic loss term. 

These can be combined via characteristic loss times B71 

giving 



L (£ p , z) =t-\E p , z) + tI(E p , z) + t~\z), 



an energy loss rate of d\nE p /dt — T^ 1 (E p ,z). We can thus 
relate the injection energy at redshift z, E' p = E p (E p , z), to 
the detected energy, E p , with 



1 dE„ 



Ep dz tt(E p ,z) dz/dt 



(8) 



The constraints imposed by this relation can be seen in Fig. 3 
of Ref. [48 1 . Using the same W(z) as above, we calculate the 
spectrum of the arriving proton flux as 



4tt 
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with dE' p {E, z)/dE p calculated numerically from Eq. (|8jl. 

In Fig. [3] we present the expected cosmic-ray proton spectra 
that result from normalizing our neutrino spectra to IceCube 
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via Eqs. (|6]i & (|7}. We note that residing near the upper edges 
of the bands would cause tension with a lack of observed PeV 
muon events. We have also checked, using CRPropa [49 1, that 
the cosmogenic gamma-ray fluxes are safely below the Fermi 
isotropic background [ 50 1 and associated constraints ll5TM55l . 

Discussion and conclusions. — If cosmic rays in the 
10 18 eV range are truly extragalactic protons, the steepness of 
the UHECR spectrum implies that their cosmic energy density 
is much larger than at the highest measured energies. If their 
escapes from acceleration regions were facilitated by photo- 
production of neutrons, then there must also be a substantial 
flux of neutrinos from pion decays. The simplified model that 
we have used captures the general flavor of the neutrino flux 
needed to explain the IceCube PeV events as a basis for com- 
parison with cosmic -ray protons. 

In Fig. [5] we see that the cosmic-ray spectrum resulting 
from our neutrino flux with ir + flavor ratios is quite close to 
measurements in the 10 18 eV range. At > 10 19 eV, the origin 
of the received flux is unclear, with a heavy-nuclear compo- 
sition inferred by Auger [25 1, in contrast to HiRes l24l . and 
the possible influence of Cen A as a local source (see ll56l ). If 
we take the Auger results at face value and assume that nuclei 
make up the difference, similar to the model of (57), a transi- 
tion in composition at ~ 10 18 5 eV results, although we have 
not addressed this or the associated neutrino fluxes in detail. 

The restrictions on this scenario include a need for accel- 
eration in the presence of a fairly soft photon background to 
result in near-threshold photoproduction. Also, a magnetic 
field of ~ 1-10 kG (depending on the T factor of the emis- 
sion) is used to cool muons, although such fields may occur 
near jet launching regions |58| or in the accretion disk, as in 
AGN core models |59 60 1. Removing the no-D requirement 
would reduce the proton flux by ~ 3; even so, it is not obvious 
that the complete set of conditions can be fulfilled. 

We also see in Fig. [3] that the cosmic -ray proton flux as- 
sociated with a neutrino flux arising from n^ 1 and fi decays 
is significantly lower than the data (and the vIp ratio could 
even be larger llSTl ). This implies that, if the requisite neu- 
trinos are made through this photopion process, the resulting 
protons are not sufficient to explain the cosmic-ray measure- 
ments. This necessitates a dominant class of accelerator that 
allows for proton escape without significant losses. 

IceCube may well be then informing us of the workings of 
a class of accelerators distinct from those providing cosmic 
rays. Dividing our inferred neutrino emissivity into sources 
yields an average of ^4— 10 x 10 43 erg s _1 (100Gpc _3 /n) 
for a space density n. These are values similar to those of 
AGN. The neutrino flux is accompanied by roughly two to 
six times this energy in gamma rays, which will quickly cas- 
cade in or near the source, and should be discussed in terms 
of concrete models (e.g., [62 -64 1. Our use of flat evolution is 
similar to that of intermediate-luminosity quasars (e.g., Il65ll ). 
Stronger source evolution, such as from bright quasars or 
the cosmic star formation rate [66 1, would result in relatively 
more neutrino flux at lower energies. This reduces the normal- 
ization and steepens the spectrum of the cosmic -ray fluxes. 



Altering the cutoffs or slopes used would not greatly change 
these basic points. Moreover, if the events turn out to be an un- 
likely fluctuation of background, the lower inferred neutrino 
flux would strengthen the conclusion that photopion interac- 
tions are not involved in freeing cosmic-ray protons from their 
sources. If not, improved data may better inform us of an as- 
pect of the extreme universe that is not evident from studies 
of cosmic rays alone. 
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